Time-reversal focusing of an expanding soliton gas in disordered replicas 
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We investigate the properties of time reversibility of a soliton gas, originating from a dispersive 
regularization of a shock wave, as it propagates in a strongly disordered environment. An original 
approach combining information measures and spin glass theory shows that time reversal focusing 
occurs for different replicas of the disorder in forward and backward propagation, provided the 
disorder varies on a length scale much shorter than the width of the soliton constituents. The 
analysis is performed by starting from a new class of reflectionless potentials, which describe the 
most general form of an expanding soliton gas of the defocusing nonlinear Schrodinger equation. 

PACS numbers: 42.65.Tg,03.75.Lm,47.35.Jk,05.10.-a 



I. INTRODUCTION 

Reversing the time dynamics of particles and waves has 
been always a central and fascinating theme in physics. 
While in classical thermodynamics irreversibility stems 
from the second principle, among several examples of 
macroscopic reversible dynamics time-reversal focusing 
(TRF) of classical waves has driven particular interest [J- 
[l3j . generalizing the concept of phase-conjugation [lij to 
non-monochromatic and transient excitation. Generally 
speaking, a TR mirror is used to detect the far-field at 
given points and time-revert the wave components prop- 
agating from a source in an inhomogeneous (scattering) 
medium, thus producing a backward propagating replica 
of the wave-packet that focuses back at the source. This 
phenomenon can occur regardless of the specific nature 
of the waves (i.e., for acoustic, electromagnetic, matter, 
elastic or seismic waves) provided dissipative effects re- 
main negligible, and has applications ranging from super- 
resolution 0, M, fici EH to enhancement of information 
transmission rate-Q, earthquake localization [l5| . atom 
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cooling [12j, image reconstruction 

Though TR have been investigated mainly in the 
framework of linear propagation, it can be extended to 
nonlinear media pi Ha. fla| . provided the governing model 
preserves time-inversion symmetry. However, in the non- 
linear case a fundamental difficulty arises from the for- 
mation of classical shock waves, which break TR locally, 
leading to irreversible loss of energy across the shock 
(where dissipation is no longer negligible). Experiments 
|16| and numerical calculations [J] show indeed that fo- 
cusing can be recovered only when TR is performed 
before the shock distance. At variance with such sce- 
nario, however, in many physical settings the regularizing 
mechanism of shocks is constituted by dispersion (rather 
than viscous effects), which leads to dispersive shock 
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waves (DSW) characterized by an expanding region filled 
with fast oscillations (modulated soliton trains) [l7l - [23l |. 
For integrable models with appropriate initial conditions 
such phenomenon can be described exactly in terms of 
a pure gas of soliton particles whose expansion becomes 
manifest beyond a point of wave-breaking determined by 
the collective behavior of solitons (20j . 

In this article we address the problem of TRF of a 
soliton gas generated beyond wave-breaking points in 
the presence of a disorder which breaks the integrabil- 
ity of the underlying model. Importantly, we explore 
the statistical stability of the TRF property, addressing 
the effect of forward and backward nonlinear propaga- 
tion occurring under different realizations of the disorder 
(with same statistical properties, i.e., two different repli- 
cas). Since a general theory of multi-solitons in weakly- 
dispersive nonlinear media subject to random perturba- 
tions is not available, nor does the methods employed 
for a single soliton give useful insight for the gas ensem- 
ble [24l - f28j . we propose a novel statistical analysis by 
combining information the ory (2i| with the replica ap- 
proach of spin-glasses [3(J El]. ^ n particular, we define 
a suitable Replica Entropy Overlap (REO) probability 
function, which acts as a thermodynamic order param- 
eter for the system. In the presence of disorder with a 
short correlation length, the REO predicts the emergence 
of a ferromagnetic phase, associated with efficient TRF 
(beyond wave-breaking points), in contrast with corre- 
lation lengths larger than the soliton width, where the 
REO yields a pseudo- glassy phase which entail spoiling 
of TRF properties. 

The paper is organized as follows. In Section [TTJ 
we introduce our model equation and analyze it in 
the limit of zero disorder, when the system is inte- 
grable. We show that a class of soliton gases can 
be characterized as reflectionless potentials in terms 
of inverse scattering transform. The shock evolution 
scenario in the integrable limit is discussed in Section 
IIII1 Section |IV] is devoted to investigate the system 
dynamics in the presence of different types of dis- 



2 



order, with both slow and fast variation (different 
spatial degrees of correlation), assessing the properties of 
TRF and contrasting it with the behavior of linear waves. 



II. A GREY-SOLITON GAS 

We focus on a universal model of dispersive wave prop- 
agation in nonlinear media, namely the defocusing Non- 
linear Scrodinger equation (NLS) |17H20| : 

^ + 0-2(H 2 -p>^2i?(,) u , (1) 

R(x) being a general integrability-breaking perturbation. 
We begin by considering Eq. (JlJ in the integrable limit 
R(x) — 0, with an input waveform 

u(x,0) = uo(x) = —p(w + iv)[wtanh(wx) + iv], (2) 

with w 2 + v 2 = 1. For p = 1, the potential uo is the well- 
known one-soliton solution of the NLS with grayness pa- 
rameter w 2 controlling the maximum depth (x = 0) and 
width of the dark (gray) notch propagating over the con- 
stant background with velocity 2v [32[ • When p ^> 1 Eq. 
J2]) contains many solitons, similarly to the limit v = 
[201 ] . However, the parameter v ^ makes the inverse 
scattering approach much more challenging than the par- 
ticular case v = 0, where the spectral transform can be 
given in closed form [20| ■ Nevertheless the essential infor- 
mation, namely the discrete spectrum (soliton content), 
can be extracted as follows. Following the approach of 
Ref. [32| we start from the linear auxiliary (Zhakarov- 
Shabat) problem of the NLS, with spectral parameter A 
(eigenvalue) : 
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= 0, (3) 



and look for its (matrix) Jost solutions ip = < / >± ( a; )i 
which are generated from the boundary values at x = 

I 



±00 obtained from Eqs. (J3j> with d/dx = and 
u Q = uo(±oo). Referring to the element <Piii x ) (the 
other elements can be found by means of symmetries 
of Eqs. © [13), and applying the coordinate change 
s = [l+tanh(u;x)]/2 as well as the F-homotopic transfor- 
mation <fy^ x = (s— l)~ l ^s~ l ~ip (we set k = \/A 2 — 4/? 2 ), 
we reduce Eq. ((3|) to the Heun equation (with poles in 
0, l,c = [w — iv)/2w, 00) for ip, defined by the following 
Riemann V symbol 33] : 



(p = V 




(4) 



where 7 = <5 = 1 — ik/2w, (3 — —ik/2w + p, a — 
~ik/2w—p, q being an accessory parameter with cumber- 
some expression 33[. Therefore, we can express <f> xl (x) 
from the full set of the 192 local solutions HI of Heun 
equation for ip. For x — ¥ —00 (s — ¥ 0), we obtain 



^ = (s iy 



l ^Hl(c,q,a,l3,j,S,s). (5) 



Then the scattering coefficients a(A) and 6(A) of the re- 
duced monodromy matrix (and hence the reflection coef- 
ficient i?(A) = b/a) should be calculated from the analyt- 
ical prolongation of to x — > +00 (s — > 1), according 
with the standard relation 132]: 



<fii 1 (s)=a(X)<p+(s)+b(\)<j>+(s). 



(G) 



However, the connection problem associated with the 
Heun equation, i.e. the analytic prolongation of Heun 
functions, lies still unsolved [33j. To overcome the prob- 
lem and arrive at the discrete spectrum, we resort to 
expanding the Heun functions HI in Eq. (|5|) in terms of 
Hypergeometric functions 2-F1 [34| : 
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E 



g n [2Fi(a,p,j-n,s)], (7) 



with g n expansion coefficients. By using the known pro- 
longation of each Hypergeometric function in s — >■ 1, we 
express the generic term in Eq. ([7]) as: 



2 -Fi(a,/3,7 - n,s) 



|r(£)l 



r(p)r(i-p) 



F + (s) + F-(s) 



rf-M 



(8) 



F ± (s) being two independent Hypergeometric functions 
near s — > 1. For a reflectionless potential containing just 
solitons both a and b go to zero for every s, a condition 
which requires the RHS of Eq. ([6]) becoming null inde- 
pendently from s. By exploiting Eq. ([8]) and well known 
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properties of Gamma functions, this condition requires: 

(9) 




with N and n positive integers. These conditions, along 
with the constraint that eigenvalues lie in the gap |A| < 
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FIG. 1. Color Online, (a) Level plot of it| 2 , showing the 
post-shock soliton expansion fans. The arrows indicate the 
two shock points as well as the slowest solitons (inner edges 
of the dispersive shock fans); (b) Snapshots of equivalent gas 
dynamic variables v and h vs. x at t = (input, dashed 
lines) and t = 0.03 (solid lines) where the left shock develops, 
the right front already exhibiting the onset of oscillations, see 
insets. The results are from the numerical integration of the 
NLS equation in the integrable limit (R — 0), with input 
uq(x), v = 0.6 and p = 40. 



2p [34], yield 2N — 1 discrete eigenvalues (solitons): 



A±„ = ±2Vp 2 -{p-nfw 2 , A = 2vp (10) 

with n — 1, N — 1 (p = N). In summary, the input 
field Uo(x) corresponds, for integer p, to a transparent 
potential describing a gas of 2p— 1 soliton particles whose 
velocities (eigenvalues) are expressed by Eq. (fTDj) . In the 
limit v = 0, Eq. (fl"0| correctly reduces to the spectrum 
already discussed in Ref. pfj| |. 



III. SHOCK SCENARIO IN THE INTEGRABLE 
REGIME 

In the strong multi-soliton regime (p = N >> 1) of in- 
terest here, the soliton particles start to become manifest 
only beyond points of wave-breaking (gradient catastro- 
phes), where sharp gradients appear in the WKB vari- 
ables: 



h(x,t) = \u\ 2 /p 2 , 

P 



v(x, t) 



(11) 



ruling the so-called dispersion- less limit of NLS [18M22] . 
These variables obey, at leading order in 1/p, hyperbolic 
equations of classical gas dynamic (h and v playing the 
role of density and velocity of the equivalent gas), which 
predict the formation of two shocks (at slight different 
times) on opposite sides of the dark notch minimum [22j . 
The integration of Eq. ([1]) confirms this scenario (see Fig. 
[T]), showing a post-shock dynamics characterized by two 
asymmetric shock fans or DSWs, which expand leftward 
and rightward, respectively. Inside these fans, solitons 
start to become manifest beyond the shock points, owing 
to regularization of the shock gradients due to dispersion. 
They appear, starting from the two slowest (and hence 
darkest) ones with eigenvalues Ao = 2Nv for the right- 
= -2y/N 2 - (N - 1) 



for v < 0) for the leftward DSW (see arrows in Fig. [TJ. 
These two asymmetric velocities define the inner edges of 
the two fans. The outer fan edges, conversely, are sym- 
metric and correspond to shallow (quasi-linear) solitons 
with large velocities A±(jv-i) = ±2\/N 2 — w 2 (these are 
not yet visible in Fig. [T|). We point out that the knowl- 
edge of the edge velocities arises, in this case, from the 
explicit solution of the direct scattering problem, while 
in non-integrable systems, in order to obtain the analo- 
gous information, one must perform lengthy calculations 
resorting to Whitham modulation theory. 



IV. TIME-REVERSAL WITH DISORDER 

Once the soliton gas propagates in the presence of dis- 
order (R 7^ 0), the expansion process of Fig. [T] is per- 
turbed with two main effects taking place: (i) the discrete 
spectrum changes (possibly including the creation of new 
eigenvalues) , as can be proven by numerically calculating 
the set Xi(t), given u(t, x) at any t (see below); (ii) a non- 
negligible amount of radiation is produced. As a result, 
the evolution could be hardly expected to be reversible 
when the wave-packet propagates backward in time over 
a different realization of the disorder. In the following 
we investigate such issue, starting in Sec. IIV Al with nu- 
merical simulations aimed at showing the TRF behavior 
of a single DSW when it propagates over a sufficiently 
long time (i.e., such that the shock fans open up consid- 
erably). After having highlighted the TRF features we 
propose in Sec. IIV Bl a statistical approach to study the 
phenomenon. Lastly, in Sec. IIV C[ we propose a simple 
approximate argument based on the Hamiltonian. 



A. Numerical results 

We have performed numerical simulations based on Eq. 
(1) with R(x) 7^ in forward and backward propagation. 
Our results show a qualitatively different behavior de- 
pending on the characteristic spatial scale (correlation) 
of the disorder R(x). A typical outcome is shown in Fig. 
[5J where we contrast by means of numerical simulations 
two qualitatively different cases: 

(i) the case of a normally distributed random variable R 
with: 




(12) 



with x being a variable coefficient, yielding correlation 
< R(x)R(x') >= a 2 S(x - x')\ 

(ii) the case of a slow (long correlation length) pseudo- 
random term which we choose as: 



ward DSW, and A_ 
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R(x) = ^g n {x- nX n ) 
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FIG. 2. Color Online. Evolution of jitj 2 for two differ- 
ent R(x) characterized by fast random (a)-(c)-(e) and slow 
pseudo-random (b)-(d)-(f) variations. The strength of disor- 
der is equal (x = 0.5) in the two cases, (a-b) color level plot 
of the forward and backward propagation; (c-d) integrabil- 
ity breaking potential R(x); (e-f) snapshots of time-reversed 
(TR) output (solid lines) compared with the input (dashed 
lines). Here the input is relative to p = 30 and w — 0.9 
(v ~ 0.44). 



with g n = a n exp[— (x/u) n ) 2 ] being Gaussian functions 
with normally distributed amplitude a n : 



< a n >= 0, 

< <4 >= <>$, 



(14) 



and displacement X n , waist uj n uniformly distributed 
in [0.15(1 — X: 0.15(1 + x))] ancl [1 — Xi 1 + xL respec- 
tively (other pseudo-random R(x) have been also tested 
and yield similar results). As shown, a remarkably good 
TRF is obtained, in spite of propagation over a different 
replica, only in the former case (delta-correlated noise), 
whereas a slowly-varying disorder leads to a poor re- 
construction with solitons scattered at random locations 
[Fig. Hb,f)]. 

The physical mechanism that leads to successful TRF 
for the case of the 5-correlated noise, is the fact that 
soliton constituents self-average the noise when the char- 
acteristic spatial scale of the latter is much shorter than 
their characteristic width. Under these conditions the 
velocities (eigenvalues) of the soliton constituents do not 
change significantly in spite of the non-perturbative char- 
acter of the random term (amplitude variations of R(x) 
are comparable with the soliton background amplitude) , 
thereby allowing to preserve the overall coherence of the 
whole soliton gas under time-reversal. This can be ap- 
preciated from Fig. [3J where we show the effect of fast 
noise [same statistical properties as in Fig. |2fb,d,f)] on 
a single dark soliton constituent, i.e. a narrow soliton 



sitting on the same background u(±oo) = p = 30 of the 
DSW modeled in Fig. O As an example we have cho- 
sen a soliton with eigenvalue Ao = 2wp, reporting in Fig. 
[3Ja,c) the evolution occurring in the presence of two dif- 
ferent realizations of the random term R(x) [displayed in 
Fig. Efb)]. From Fig. [3jc) it is clear that the net ef- 
fect of the random term is to add small random fluctua- 
tions to the soliton background and a small displacement 
in the return position. Such features, being due to the 
statistical difference of the fluctuations in forward and 
backward propagation, disappear when backward prop- 
agation occurs under the same statistical realization of 
the noise, which leads to perfect reconstruction as shown 
for comparison in Fig. [3jd). However, even in the pres- 
ence of different statistical realizations, we observe that 
the net change of velocity of the soliton (which has ob- 
viously statistical nature) remains extremely small, e.g. 
in Fig. |3jb) we measure a net change Av/v < 0.4%. 
This allows for recombining all the solitons at the right 
locations thus reproducing the reversal of the entire gas 
or DSW shown in Fig. [2a). It is important, however, 
to emphasize that it is only the combination of two in- 
gredients, namely fast fluctuations and soliton features, 
that makes possible to obtain a rather good TRF. In gen- 
eral non-soliton wave-packets do not share this feature. 
For instance, in the linear regime, where the evolution is 
dominated by linear spreading (diffraction), TRF is not 
observed. In order to show this, we report in Fig. [4] the 
outcome of a simulation of the evolution of a linear wave- 
packet. For comparison we inject the same waveform as 
in Fig. [3J though with a hundred times lower |m| 2 , and 
evolve the wave-packet in the same medium (i.e., same 
statistical properties of R(x)), over the same temporal 
scale. As shown in Fig. [4ja,c), time- reversal turns out 
to be quite poor, with fluctuations which become compa- 
rable with the main dark notch injected in the medium. 
This must be contrasted, once again, with perfect time- 
reversal reconstruction which is observed when backward 
propagation occurs under the same realization of noise, 
as reported in Fig. |4jb,d) for comparison. 

Conversely, when the characteristic length scale of 
R(x) becomes comparable or larger than the width of 
soliton constituents of the DSW, scattering of the solitons 
occurs according to the local features of the noisy poten- 
tial R(x). In this regime, significant statistical changes of 
the soliton velocities become important and lead to ham- 
per the TFR feature. In order to deepen this point we 
have calculated the distribution of soliton velocities by re- 
sorting to numerical calculation of the eigenvalues of the 
spectral problem [Eqs. (J3J)], starting from the potential 
u(x, t) obtained numerically. To this end we have em- 
ployed an algorithm based on Rayleigh-Ritz variational 
procedure [351 ] . that we have adapted to deal with the 
defocusing case (self-adjoint scattering problem). In Fig. 
[5jb) we report the eigenvalues A„ (see triangles) vs. a 
progressive integer n that identifies soliton components 
(soliton number) of the total field [shown in Fig. [Ha)] 
emerging after forward propagation (at t=0.1), in the 
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FIG. 3. Color Online. Evolution of a 1-soliton with back- 
ground p = 30, and eigenvalue Ao = 2vp (w = 0.9, v — 
vl — to 2 ), in the presence of a fast random term R(x) [same 
statistical properties as in Fig. [2] (left column)], (a) color 
level plot of the forward and backward propagation; (b) inte- 
grability breaking potential R(x); (c) snapshots of TR (solid 
red line) and forward (thin solid black line) output, compared 
with the input (dashed blue line). The inset shows the net 
displacement occurring upon TRF; (d) reversed waveform for 
equal R(x) in forward and backward propagation (ideal phase- 
conjugation): the TR output is exactly superimposed to the 
input. 




FIG. 4. Color Online. Evolution of \u\ for a linear wave- 
packet (same waveform as in Fig. 3, 100 times weakened 
\u\ 2 ) under 5-correlated R{x) [same statistical properties as 
in Fig. 2(left column) and Fig. 3]. (a-c) TRF when back- 
ward propagation takes place under a different realization of 
noise [R(x) not shown, similar to those of Fig. [3jb)]; (b-d) 
ideal phase conjugation for the same realization R(x) in for- 
ward/backward propagation. 
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FIG. 5. Color Online. Case of slow pseudo-random R(x) 
[same statistical properties as in Fig. Elb-d-f)]: ( a ) snapshots 
of the rightward expanding fan after forward propagation at 
t = 0.1, comparing the noisy case (red solid curve) with the 
integrable case (black solid curve, R(x) = p). The dashed blue 
curve stands for the input, (b) relative soliton eigenvalues X n 
(on the abscissa) vs. soliton number n in the right half of the 
gap, contrasting the integrable case (open circles) with the 
noisy one (triangles). 



presence of R(x) with long correlation length (as in Fig. 
HI right column). These data are contrasted with the 
situation such that random variations are switched off 
(open circles in Fig. [U[b); for sake of comparison with 
the noisy case, we maintain R(x) equal to the mean value 
p = 30 and use the same algorithm, though, obviously, 
the outcome is the same as in our analytical scattering 
approach in the integrable case, due to iso-spectrality 
and phase-rotation symmetry of the NLS equation). By 
comparing the two sets of data reported in Fig. EJa,b), it 
is clear that significant random changes of velocities take 
place, especially for what concerns darker solitons, which 
are expected to contribute more to the reconstruction of 
the input notch profile. Also several new eigenvalues are 
created primarily at the edge of the spectral gap, where 
solitons turn out to be denser. Such alteration of the dis- 
tribution of soliton velocities (and width and darkness as 
well, since they are related to velocity) is the clear signa- 
ture of the fact that the noise is not simply averaged out 
in this case, thus spoiling TFR. 



B. Replica theory: 



ferromagnetic vs glassy soliton 
phase 



To proceed further and characterize the process of TRF 
statistically, we provide an approach based on the replica 
theory of spin-glasses. We consider two different observ- 
ables of the system, {/*•"'}, {/^}, emerging from two 
different TRF experiments labeled a and /3. Being inter- 
ested in the TRF reproducibility of the dark notch, we 
pose / = \p — | it 1 1 (we remove the background p), nor- 
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malizing so that f^°° dxf(t,x) = 1. Then we define a 
suitable entropy distance q a ^ between {f^ a >}, {f^}' 

[q*,?]- 1 = 1 + \[K{a, (3) + K(0, a)}, (15) 

K(a, (3) = J dxf^ log jf^j being the (non-negative) 
Kullbach information (29[. By construction, < q a p < 
1, with the extremes attained when the states a and 
/3 turn out to be identical (q a ,p = 1 or ferromagnetic 
phase), or viceversa uncorrelated (q a 8 = 0, equivalent 
to a paramagnetic phase). Then in order to describe the 
statistics of all possible distances in the space of states, 
we introduce the probability distribution function P{q, t): 

P{q,t)=< S(q a ,p-q)>, (16) 

with < ... > denoting an ensemble average. It is the 
shape of this function P(q,t), which we call replica en- 
tropy overlap, that plays the role of a physical order pa- 
rameter for the system. Figure [6] shows the behavior of 
P(q,t), as calculated by averaging over 20 samples with 
different degree of disorder (i.e. \ € [0.01,0.5]), each 
characterized by 20 different specific realizations of R(x). 
The number of samples and quenched replicas have been 
chosen in order to let P(q, t) converge to a steady-state 
distribution. As shown, for a noise with short correla- 
tion length, the evolution exhibits a tendency to remain 
in a ferromagnetic phase characterized by a probability 
function peaked around q ~ 1. In this case the soli- 
ton wave-functions "sample" an average disorder during 
propagation and are therefore able to revert the evolu- 
tion even in the presence of a different replica. Viceversa, 
when the disorder R has long correlation length, the tem- 
poral evolution of the system exhibits a clear tendency 
towards what is known as a glassy phase |3lj . the latter 
being characterized by a space-of states filled continu- 
osy with all the possible distances in [0, 1]. In this case, 
irreversibility becomes manifest by the fact that, upon 
backpropagation, P(q, t) no longer return to be peaked 



Equations (|T7| - (|18|) entail that the dynamic is that of an 
ensemble of kinetic free particles, perturbed, at leading 
order, by a term which depends only on the spatial mean 
of the disorder potential. The motion ruled by Eqs. (fl7|) - 
(|18p is completely reversible and, more important, does 
not depend on the particular realization of R. This re- 
sult relies on nature of the soliton wave- functions that, 
being more extended than the characteristic spatial scale 



at q = 1, but continuosly fills the interval Aq g [0.5, 1] 
(Fig.©. 



C. Hamiltonian perturbative analysis 

In order to gain more physical insight on such be- 
havior, we provide a simple though approximate ar- 
gument. Starting from the Hamiltonian density % = 
\u x \ 2 + (\u\ 2 - p 2 ) 2 + 2R\u\ 2 of Eq. flT), and assuming 
the randomness to be fast varying R = R(x/e) (e <C 1) 
and such that the form of each soliton is not appreciably 
modified, the Hamiltonian H = J %dx becomes: 




0.5 1 0.5 

q q 



FIG. 6. Color Online. Countour plot of entropy overlap 
probability P{q) vs. time, obtained with disorder R(x): (a) 
uncorrelated (fast) random term with normal distribution; (b) 
correlated pseudo-random term obtained through summation 
of Gaussians (see text); The probability has been calculated 
by mediating over 20 samples (x £ [0.01, 0.5]), each simulated 
with 20 different realizations of disorder. The input has been 
taken as a dense soliton gas with p = 30 and w = 0.9. 

where the sum represents a purely kinetic term given 
through the unperturbed Hamiltonian of the single soli- 
ton Hi = | (Ap 2 — A|) 3 ^ 2 . The evaluation of the last 
integral in Eq. (|17[) can be made in the spatial frequency 
domain as follows (hat denoting Fourier transform): 



(18) 

I 

of R, do not see the landscape of R but just its average. 
This explains in simple terms why TRF occurs regard- 
less of the specific configuration of disorder. In the case of 
slowly varying terms, conversely, Eq. (|18[) no longer holds 
and the soliton wave-functions evolve according with the 
valley and peak structure of R, thus spoiling statistical 
stability of TFR against different replicas. 



lim 



\u\ 2 R 



-L 



e 1 ' x dxdk = 2ep lim 
e / L— >-oo 



Ft hT 

ft S i n ^dk + 0{e 2 ) = 2ep 2 R(0) + 0(e 2 ). 
k e 
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V. CONCLUSIONS 



In summary, we have presented an example of a 
reversible shock dynamics by considering the refocusing 
of an expanding soliton gas propagating over different 
replicas of a disordered channel. A successful recompres- 
sion of the gas occurs for a fast disorder, thanks to the 
onset of a ferromagnetic phase in the order-parameter 
which characterize the process. The most obvious exper- 
imental test or our predictions can be made in nonlinear 
optics where DSW have been observed (though in 



settings better described by generalized NLS equations 
21]) and R(x) represents a random refractive index 
distribution. However, the universality of the NLS model 
let us envisage also experimental tests in other fields, for 
instance with matter waves in Bose-Einstein condensates. 
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